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LONG-TERM  GOALS 

The  long  term  goals  of  this  project  are  the  development  and  application  of  new  methods  of 
investigation  for  the  use  of  Lagrangian  data.  Special  attention  is  given  to  the  combined  use  of  data 
and  models,  with  focus  on  applications  to  regional  and  coastal  areas. 

OBJECTIVES 

1)  To  develop  and  apply  new  techniques  for  the  assimilation  of  Lagrangian  data  in  ocean  gen¬ 
eral  circulation  models  (OGCMs). 

2)  To  improve  previous  results  on  statistical  studies  of  particle  motion  using  data  and  stochastic 
models. 

APPROACH 

The  work  involves  a  combination  of  numerical,  analytical  and  statistical  techniques. 

WORK  COMPLETED 

1)  Development  of  a  Lagrangian  data  assimilation  method  and  its  application  to  a  quasi-geostrophic 
model.  The  results  are  accepted  for  publication  in  J.  Geophys.  Res. 

2)  Analysis  of  drifter  data  and  numerical  results  with  application  to  large  scale  ocean  flows 
(Tropical  Pacific)  and  regional  coastal  flows  (Sicily  Channel).  Results  are  published  in  J. 
Geophys.  Res.  and  J.  Mar.  Sys. 

RESULTS 

Work  in  the  present  grant  period  has  been  focused  on  two  main  topics: 

1)  A  specific  formulation  for  assimilation  of  Lagrangian  data,  which  takes  directly  into  ac¬ 
count  the  nature  of  the  measurements,  i.e.  the  positions  of  drifting  buoys,  has  been  tested  in  the 
framework  of  a  double-gure  quasi-geostrophic  (qg)  model  and  of  a  primitive  equation  (MICOM) 
model 

A  central  problem  in  the  assimilation  of  Lagrangian  data  is  that  there  is  a  nonlinear  relationship 
between  the  observed  float  positions  and  the  model  variables  to  be  modified,  i.e.  the  Eulerian 
velocities.  In  order  to  address  this  issue,  a  new  formulation,  in  the  framework  of  the  Optimal 
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Interpolation  (01)  method,  has  been  implemented  in  collaboration  with  L.  Piterbarg  (Molcard  et 
al.,  2002).  Conceptually,  the  drifter  position  is  forecasted  by  the  model  during  the  assimilation 
interval,  and  the  model  Eulerian  velocity  is  corrected  in  order  to  minimize  the  distance  betwen 
forecasted  and  observed  positions. 

A  simplified  version  of  the  formulation  (local  in  time  in  space)  has  been  tested  in  the  framework 
of  a  double-gyre  qg  model.  An  extensive  study  has  been  conducted  to  quantify  the  effectiveness 
of  the  assimilation  using  the  twin  experiment  approach.  The  method  is  found  effective  in  a  vaste 
range  of  parameters,  in  terms  of  number  of  drifters,  frequency  of  assimilation  and  uncertainties 
associated  with  the  forcing  functions  driving  the  ocean  model.  The  performance  of  the  Lagrangian 
assimilation  technique  is  also  compared  to  that  of  conventional  methods  (“pseudo-Lagrangian”)  of 
assimilating  drifters  as  moving  current  meters  (e.g.,  Ishikawa  et  al.,  1996).  The  fully  Lagrangian 
assimilation  provides  better  results  in  all  the  tested  cases  (Fig.  1). 

The  performance  of  the  assimilation  technique  is  presently  tested  in  the  framework  of  prim¬ 
itive  equation  models,  using  MICOM.  A  simultaneous  correction  of  the  layer  thickness  is  also 
performed,  using  a  dynamical  relationship  based  on  geostrophy.  Results  are  very  encouraging  as 
shown  in  the  example  of  Fig.  2. 


(a)  (b) 
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Figure  1:  (a)  Typical  convergence  error  using  fully  Lagrangian  (solid  line)  and  pseudo-Lagrangian 
(dashed  line)  assimilation  schemes,  (b)  Normalized  error  difference  between  Lagrangian  and 
pseudo-Lagrangian  methods  as  a  function  of  assimilation  period.  The  normalized  error  differ¬ 
ence  between  the  fully  Lagrangian  method  and  assimilation  of  current  meter  data  is  marked  by 

(  ❖  )  .  The  Lagrangian  time  scale  is  10  days. 
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2)  Study  of  transport  properties  and  statistics  of  particle  motion  in  the  Adriatic  Sea  using  a 
drifter  data  set  ( 1990-1999). 

The  work,  performed  in  collaboration  with  P.  Poulain  and  A.  Maurizi,  builds  on  previous  re¬ 
sults  (Falco  et  al.,  2000)  obtained  using  a  more  restricted  data  set  (1994-1996).  The  use  of  the 
extended  data  set  (Poulain,  2001)  allows  to  explore  the  effects  of  inhomogeneity  and  nonstationar- 
ity  (Maurizi  et  al.,  2002)  which  were  neglected  in  Falco  et  al  (2000).  The  Adriatic  Sea  have  been 
partitioned  into  3  zones  based  on  eddy  kinetic  energy  and  topography,  and  the  properties  of  the 
turbulent  velocity  have  been  studied  separetely  in  each  zone.  In  Fig.  3,  the  velocity  autocorrelation 
of  is  shown.  As  it  can  be  seen,  at  initial  times,  the  autocorrelations  are  similar  in  the  3  zones  and 
approximately  exponential,  while  at  longer  times  the  autocorrelation  in  zone  2  shows  a  signifi¬ 
cantly  different  behaviour,  and  it  is  characterized  by  a  longer  time  scale  of  the  order  of  15  days, 
reflecting  also  on  the  global  average.  Zone  2  corresponds  to  the  most  energetic,  boundary  coastal 
area  in  the  southern-central  Adriatic,  which,  as  shown  by  Falco  et  al.  (2000),  is  the  most  difficult 
zone  to  model  in  terms  of  transport  properties.  The  existence  of  the  secondary  scale  explains  these 
difficulties  and  suggests  a  more  complete  transport  parameterization. 


Figure  3.  Log  of  turbulent  velocity  autocorrelations  R(t )  computed  in  the  3  zones  of  the  Adriatic 
Sea  and  over  the  whole  basin  (global).  The  three  steepest  lines  indicates  fits  to  data  in  the  3  zones 
at  initial  times  ( approximately  first  3  days),  while  the  fourth  line  indicates  fit  to  data  at  longer  times 
for  zone  2.  The  slopes  indicate  e-folding  time  scales  x.  All  zones  have  initial  time  scales  of  order 
of  1.5  days.  In  zone  2,  there  is  a  second  time  scale  of  15  days  for  global  average. 

IMPACT/APPLICATIONS 

The  assimilation  of  Lagrangian  data  is  expected  to  have  a  direct  impact  on  ocean  prediction 
using  OGCMs.  The  turbulent  transport  studies  provide  tools  for  parameterization  and  transport 
modeling. 

TRANSITIONS 

The  Lagrangian  assimilation  results  are  planned  to  be  used  in  OGCMs,  such  as  a)  MICOM/HYCOM, 
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in  collaboration  with  L.  Piterbarg,  A.  Mariano  and  M.  Chin;  and  b)  OPA  model  of  the  Mediter¬ 
ranean  Sea  in  collaboration  with  N.  Pinardi.  The  transport  study  in  the  Mediterranean  Sea  will  be 
carried  out  in  collaboration  with  P  Poulain  and  E.  Zambianchi. 

RELATED  PROJECTS 

(i)  Statistical  problems  in  ocean  modelling  and  prediction.  PI:  L.  Piterbarg,  ONR 

(ii)  Mediterranean  forecasting  project,  PI:  N.  Pinardi,  EU 
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